The Thermopower of Quantum Chaos 
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The thermovoltage of a chaotic quantum dot is measured using a current heating technique. 
The fluctuations in the thermopower as a function of magnetic field and dot shape display a non- 
Gaussian distribution, in agreement with simulations using Random Matrix Theory. We observe no 
contributions from weak localization or short trajectories in the thermopower. 
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The electrical conductance of small - characteristic size 
much smaller than the electron mean free path - con- 
fined electron systems (usually denoted as quantum dots) 
shows distinct fluctuations. These fluctuations display 
correlations as a function of an external parameter such 
as shape or magnetic field, which can be described in a 
statistical manner. The electrons can, in fact, be viewed 
as billiard balls moving in a classically chaotic system 
where many random reflections at the system walls oc- 
cur. Because of the wave-like nature of the electrons, 
quantum mechanics is needed to describe these systems 
fully. Chaos in quantum dots has been investigated 
in conductance measurements but the analysis turns out 
to be difficult. So-called short trajectories Q and weak 
localization effects add up to the signature of chaotic 
motion. Moreover, current heating of the electrons in the 
dot appears to be unavoidable in conductance measure- 
ments. Electron heating effects in the dot smear out the 
underlying chaotic statistics and therefore the observed 
fluctuations exhibit mostly a Gaussian distribution, al- 
though theory predicts non-Gaussian distributions when 
a small number of electron modes is admitted to the dot 
H . Only when dephasing (modelled as extra modes cou- 
pling the dot to the environment) is included, Random 
Matrix Theory (RMT) gives a Gaussian distribu- 
tion. Very recently, Huibers et al. || observed small 
deviations from a Gaussian distribution in conductance 
measurements. However, other transport properties cal- 
culated from these data exhibit again Gaussian distribu- 
tions in contrast to theoretical predictions. 

An alternative for the conductance measurements pur- 
sued so far (which inherently are accompanied by electron 
heating inside the dot) is to investigate the thermoelec- 
tric properties of a system. Thermopower measurements 
have already been used to study semiconductor nanos- 
tructures like quantum point-contacts J9] and quantum 
dots in the Coulomb blockade regime m|,0|. The ther- 
mopower S measures directly the parametric derivative 
of the conductance, S oc G~ 1 dG/dX with X = E (en- 
ergy), and thus yields both similar and additional in- 
formation on the electron transport processes as can be 



obtained from conductance measurements. The distri- 
bution of parametric derivatives (X — E,B, shape, . . .) 
of the conductance of a quantum dot is the subject of 
recent RMT- investigations fl2| , [l3|| . The probability dis- 
tribution for the thermopower is again expected to be 
non-Gaussian for chaotic conductors, exhibiting cusps at 
zero amplitude and non-exponential tails [ p^| , pd[ . 

In this paper, we present magneto-thermopower mea- 
surements of a statistical ensemble of chaotic quantum 
dots. The observed thermopower fluctuations show a 
non-Gaussian distribution. We present a numerical fit 
based on RMT which describes the experimental data. 
We demonstrate that effects like short trajectories, weak 
localization and dephasing are absent in thermopower 
measurements . 

In Fig. |l|a the measured device is shown schematically. 
A quantum dot (lithographic size 800 nm x 700 nm) is 
electrostatically defined (gates A, B, C and D) in a stan- 
dard, high-mobility 2-dimensional electron-gas (2DEG) 
in a GaAs-(Al,Ga)As heterostructure. The 2DEG has a 
mobility of /x ~ 10 6 cm 2 (Vs) _1 for an electron density 
of 3.4 x 10 11 cm~ 2 at 4.2 K. A 2 fim wide and a 20 jum 
long electron-heating channel is defined next to the quan- 
tum dot (gates A, D, E and F). The sample is kept at 40 
mK in a dilution refrigerator equipped with a supercon- 
ducting magnet. Transport measurements are performed 
using standard phase-sensitive techniques. For reasons of 
comparison all data shown in this paper were obtained 
from the same sample. We have obtained similar results 
in several other devices. 

Conductance data are shown in Fig. |]b. The graph is 
the magnetoresistance of the dot, averaged over a large 
number of different configurations |lq] . Both point con- 
tacts leading to the dot were adjusted to a conductance 
G — 4e 2 /h corresponding to two spin-degenerate modes 
in the point contacts. As in Ref. ||, an ensemble of 
configurations was created by repeatedly changing the 
voltage on gate B by a small amount (8V^ = 10 mV). 
As is evident from the figure, apart from chaotic conduc- 
tance fluctuations also the signatures of weak localization 
(sharp peak around B = 0T) and short trajectories 
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FIG. 1. (a) Schematic topview of the measured sam- 
ple. The crosses denote the Ohmic contacts to the 2DEG; 
the hatched areas denote the gates. The heating cur- 
rent is applied between II and 12. The thermovoltage is 
measured between VI and V2. The quantum dot is de- 
fined by applying a negative voltage to gates A,B,C and 
D. (b) Magnetoconductance of the dot, averaged over a 
large number of different dot configurations. The influ- 
ence of short trajectories is characterized by the dashed 
line. Inset: Conductance distribution for \B\ > 50 mT. 



(characterized by conductance background that ex- 
hibits a polynominal dependence on magnetic field [pB . 
Fig. lb, dashed line) are clearly observed. In order to ex- 
tract the statistics of the fluctuations, the conductance 
measurements were corrected for these features. The 
resulting distribution of the conductance fluctuations is 
shown in the inset of Fig. [l]b (for magnetic fields larger 
than 50 mT) . A Gaussian function fits these data well. 

By passing a low-frequency (13 Hz) current through 
the ohmic contacts Ii and I2, the electron gas in the chan- 
nel is heated (T e oc I 2 ) while the wide 2DEG regions re- 
main in equilibrium with the lattice. The ac-heating cur- 
rent / w 0.4 //A is small enough to avoid lattice heating 
and we have ensured that we are in the regime of linear 
response. The temperature difference between heating 
channel and electron reservoir induces a thermovoltage 
Vth across the chaotic dot. Vth can be measured, using 
quantum point-contact (QPC) E-F as a reference point, 
between ohmic contacts Vi and V2. Vth is then related 
to the thermopower of the dot, 5dot, as 

V th = V 2 -V 1 = Vd ot - V re f = (Sdat - Srcf)AT (1) 

Here, QPC E-F is adjusted such that its thermopower, 
SVof, is minimal and constant for all measurements. As 
in the conductance measurements, the transmittance of 
QPCs A-D and C-D was adjusted to G — 4e 2 /h. Again, 
varying the voltage applied to gate B, V^, was used to 
change the shape of the dot. 



In Fig. ||a a grayscaleplot is shown of the transverse 
voltage Vth as a function of magnetic field and V® . The 
step in gate voltage between two successive magnetic- 
field sweeps was 8V® = 10 mV. The magnetic field range 
was limited to \B\ < 150 mT to avoid the regime where 
the quantum Hall effect becomes dominant. The char- 
acteristic fluctuations are stable in time and well repro- 
ducible. As an example, the trace for Vj 8 = —550 mV is 
plotted separately in Fig. 2b. The fluctuations are sym- 
metric around B = 0T with a zero mean. 

Because the resistance of the electron heating chan- 
nel is magnetic-field dependent due to classical (break- 
down of the entrance resistance) and quantum (weak 
localization and conductance fluctuations) transport ef- 
fects, also the electron temperature in Eq. 1 is some- 
what dependent on magnetic field, T e = T e (B). For 
the given device structure, this dependence can easily 
be determined experimentally, using the quantized ther- 
mopower of a point contact J^]. The thermovoltage 
of QPC A-D, adjusted for maximum thermopower, is 
measured as a function of magnetic field while keep- 
ing V® = Vg — V, i.e. without defining the dot, 
yielding Vth .channel = (Sad — S Ic i)AT(B). The varia- 
tion in channel-temperature, which turns out to be only 
a few percent, is effective for each individual measure- 
ment. Thus, by dividing the dot thermovoltage by the 
QPC thermovoltage, this effect can be eliminated. We 
have verified that this calibration does not influence the 
statistics of the thermopower fluctuations. 

The inset of Fig. |^b shows the temperature-variation 
corrected thermovoltage, averaged over all configura- 
tions. We show this averaged trace to illustrate that in 
contrast to the averaged conductance (Fig. lb) signatures 
of weak localization and short trajectories are absent in 
the averaged thermopower jlj|. For the weak localiza- 
tion correction, this is readily understood when one re- 
alizes that the correction to the conductance for zero di- 
mensional system is energy independent and any signs of 
weak localization in the thermopower should derive from 
the energy dependence of the phase-coherence length 1$ 
which is presumably small. (However, see Ref. [|l6| for ex- 
periments on weak localization thermopower in the two- 
dimensional quantum-diffusive transport regime). The 
absence of a signature of short trajectory effects in the 
thermovoltage measurements implies also a very weak 
energy dependence for the conductance correction due to 
these processes. Heuristically one might argue that the 
fast transit times involved with short trajectories cor- 
responds to large energy scales while the thermopower 
measures a local (at Ep) derivative. It would be inter- 
esting to investigate these effects theoretically. 

We now proceed to compare the statistics of the ob- 
served thermopower fluctuations with theoretical predic- 
tions. The system symmetry, denoted in RMT with an 
integer /3, changes with magnetic field; around B = T, 
time-reversal symmetry (TRS, (3 = 1) is present while for 
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FIG. 2. (a) Grayscaleplot of the thermovoltage of the 
quantum dot as a function of magnetic field (X- Axis) and 
of gate voltage on gate B (Y-Axis). The gate voltage is 
changed by a constant small amount of 5V^ = 10 mV for 
each magnetic field sweep. The light areas denote a large 
(maximum 5.5 fiV) positive thermovoltage, the dark ar- 
eas a large (maximum -5.5 fiV) negative thermovoltage. 
(b) Individual thermovoltage trace for = —550 mV 
[dashed line Fig. (a)]; Inset: Thermovoltage, corrected 
for the temperature variation of the heating channel and 
averaged over a large number of dot configurations. Sig- 
natures of weak localization and short trajectories are 
absent. 



higher magnetic fields this symmetry is broken {(3 = 2). 
The transition between the two regimes is a gradual one; 
for the present analysis we restrict ourselves to the two 
extremes. Experimentally, the regime with TRS is the 
magnetic field range where the weak localization effect 
dominates conductance measurements (\B\ < 30 mT), 
TRS is broken for \B\ > 50 mT. Counting the fluctua- 
tion amplitudes of the corrected thermopower for these 
two regimes leads to the histograms presented in Fig. |^. 
The dashed lines are the best Gaussian fit to the exper- 
imental data. Clearly, strong deviations from Gaussian 
statistics occur. 

These deviations are also expected from RMT calcu- 
lation. In Ref. |L4| the thermopower distribution of a 
chaotic quantum dot has been obtained for single-mode 
contacts. The distribution exhibits a cusp at S = and 
tails as P(S) oc jS 1 ! -1- ^ In \S\, which displays a clear de- 
viation from a Gaussian distribution. However, for a 
large number of conducting channels this distribution be- 
comes again Gaussian. For the experimental data, taken 
with two spin-degenerate conducting modes in the leads, 
deviations from the a Gaussian distribution are still ex- 
pected. Since in RMT an analytical treatment for the 
thermopower distribution is possible only for single-mode 
leads, Monte-Carlo simulations have to be employed for 
the present system. 
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FIG. 3. (a) Thermopower distribution in presence of 
TRS < 40 mT). Experimental results (dots), simu- 
lation results (solid line) and Gaussian fit (dashed line) 
(b) Thermopower distribution for broken TRS (\B\ > 50 
mT). Experimental results (dots), simulation results 
(solid line) and Gaussian fit (dashed line). Inset: Experi- 
mental thermopower distribution (dots) for broken TRS, 
compared to simulation results including strong dephas- 
ing (solid line). 



The Hamiltonian 7i of the closed dot is drawn ran- 
domly from the Gaussian ensemble 



P(H) oc e 



-cTrWH 



(2) 



where c is a constant setting the mean level spacing A 
and Ti is real symmetric (/3 = 1) or hermitian ((3 = 2). 
The scattering matrix fl% M of the open system is cal- 
culated from 7i using |18|[ 

M(E) = 1 - 2mW i (E-H + iixWW^) ~ l W. (3) 

Here, W is a rectangular matrix coupling the states in 
the dot to the scattering channels. The thermopower is 
then calculated from 



5' 



r 2 klT d 



e dE v ' 



(4) 



E — E~p 



where T(E) = £ Q 



is the probabil- 



ity for the transmission from lead 1 to lead 2. The 
differentiation is done numerically for each realization 
of the Hamiltonian. The density of states p(Ep) = 
l/2ni Tr (M^ dM / dE)\ E=EF is used as a weight factor to 
account for a large charging energy E§. The resulting 
distributions of the thermopower fluctuations for both 
symmetry classes are shown as solid lines in Fig. |[ Evi- 
dently, the simulations represent the experimental results 
much better than a Gaussian distribution function. 
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For a comparision of experimental and theoretical dis- 
tributions, the horizontal axes of both data sets have to 
be scaled, while keeping the normalization (area equals 
one). Only one scaling parameter was needed to fit the 
distributions for TRS and broken TRS (Fig. 3), which im- 
plies a temperature difference across the dot of AT = 56 
mK. As an independent check for the correctness of this 
scaling procedure, we have in addition measured the ther- 
movoltage (for the same heating current) across QPC 
AD, adjusted for maximum thermovoltage, i.e. between 
the N = 1 and N — 2 conductance plateaus ||. For this 
configuration, the thermopower of a QPC is quantized, 
directly yielding a value for AT. We find AT = 54 mK 
at I = 0.4 fiA, in good agreement with the scaling result. 

It is possible to show that in thermopower measure- 
ments dephasing corrections, which are used to explain 
Gaussian conductance fluctuation distributions of chaotic 
quantum dots, are indeed irrelevant. The influence of a 
finite phase-coherence time can be studied by con- 
necting the dot to a third virtual reservoir. The bath is 
coupled via N$ 1 channels, corresponding to the in- 
coherent limit (h/r^ 3> A). The chemical potential /i^ 
and temperature of the bath are chosen exactly in 
between those of the real reservoirs, such that it draws 
neither charge nor heat current. The distribution of the 
thermopower for broken TRS (/? = 2) then given by 

P(a) = 2--^A^M 4 ff ( 8 ^ 2 7 fc ) ±\2Na 

k=0 ^ ' ' 

The result of Eq. is plotted for N = 2 in the inset 
of Fig. ||b (solid line)! It can be seen that there is only 
little agreement between the experimental thermopower 
data and theory including dephasing which proofs that 
the observed statistics do not indicate the presence of any 
dephasing induced by electron heating. 

To conclude, we have demonstrated that the ther- 
mopower measurements on a chaotic quantum dot re- 
veal directly the theoretically predicted non-Gaussian 
fluctuation distributions. In contrast to the best re- 
sults obtained so far in conductance measurements, we 
don't have to include thermal broadening or dephasing 
to model our experimental results. Another distinction 
from conductance measurements is that the thermopower 
data are not influenced by weak localization or short tra- 
jectories. Therefore, thermopower measurements can be 
considered as an excellent tool in the area of investigating 
chaotic quantum transport properties in open systems. 
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